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ABSTRACT:  This  report  presents  the  results  of  ballistics  range  and 
wind-tunnel  measurements  -of  the  static  force  and.  moment  and  the  pitch¬ 
damping,  roll-damping  and  Magnus  moments  for  a  General  Purpose 
20-Millimeter  Projectile.  Data  are  presented  for  a  Smooth  and  Basic 
configuration.  These'  aerodynamic  data  are  used  in  a  special  digital 
computer  program  to  rapidly  assess  the  stability  of  the  projectiles 
over  various  mission  profiles.  The:  normal-force  and  pitching-moment 
derivatives  are  evaluated,  theoretically,  by  Wood's  method  and 
compared  with  the  measurements. 
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INTRODUCTION: 

The  purpose  of  this  report  is  twofold:  first,  aerodynamic  static 
and  dynamic  measurements  are  presented  for  the  Basic  20mm  general 
purpose  projectile  and  for  configurations  -which  are  simple  modifications 
of:  this  basic  shape;  secondly,  these  aerodynamic  data  are  used  in 
a  , specially  developed  computer  program  for  the  rapid  evaluation  of 
projectile  stability.  Details  of  this  program  are  presented  in  this 
report. 

This  investigation  program  was  carried  out  in  several  phases i 
Initially,  a  number  of  configurations  were  selected  from  simple 
modifications  of  the  basic  shape.  A  simple  modification  is  defined 
as  a  single  change  in  an  important  feature  of  the  shape,  such  as 
smoothness,  length,  base  geometry,  etc.  The  change  of  two  or  more 
of  these  salient  configuration  features  would  be  considered  a 
compound  nusdifigat ion.  In  this  early  phase  these  modifications 
consisted  of  a  0.6-  and  0.3-caliber  length  extension,  a  seven- 
degree  0.5-baliber-lcng  boattail  and,  finally  a  smooth  body.  This 
last  shape  was  formed  from  the  Basic  body  by  omitting  the  rotating 
band’  and  filling  in  the  crimping  and  fuze  grooves.  The.  Basic 
configuration,  which  will  be  discussed  in  more  detail  in  the 
next  section,  is  a  4 . 9-caliberr-long:  20mm  projectile  intended  for 
general  use  within  the  U.  S.  Navy. 

As  a  result  of  this  initial  investigation  it  was  decided,  for 
various  reasons,  to  ^retain  the  Basic  configuration  as  originally 
formulated.  However,,  it  was  felt  that  it  would  be  of  value  to 
obtain  a  fairly  ^detailed  set  of  aerodynamic  data  on  the  smooth 
body.  If  removable,  or  consumable,  rotating  bands  have  wide  use 
and  if  improved  projectile  case  bonding  eliminates  the  crimping 
groove,  the  low-drag  advantages  of  the  smooth  configuration  might 
be  approached. 

These  tests  ’/ere  carried  out  in  the  Naval  Ordnance  Laboratory's 
Supersonic  Tunnel  No.  1  and  Pressurized  Ballistics  Range.  The 
wind  tunnel  was  used  to  obtain  nonlinear  variations  of  the  normal- 
force  and  pitching-mpment  coefficients  and  the  Magnus  force  and  moment 
coefficients  with  angle  of  attack.  The  ballistics  range  tests  were 
intended  to  obtain  the  damping-in-roll  and  the  damping-in-pitch 
derivatives  which  were  not  available  from  the  wind-tunnel  measurements. 
The  ballistics  range  also  furnished  the  normal-force,  pitching- 
moment  and  Magnus  derivatives,  providing  a  common  ground  for 
comparison  between  the  wind  tunnel  and  ballistics  range. 

In  addition  to  the  aerodynamic  data,  this  technical  report 
presents  a  digital  computer  program  assessment  of  the  stability  of 
the  configuration  under  a  variety  of  conditions.  The  program  is 
presented  in  detail  in  one  of  the  sections.  Finally,  a  comparison 
is  made  between  the  wind-tunnel-ballistics  range  measurements  of 
the  normal-force  and  pitching -moment  derivative  with  an  analytic 
evaluation  of  these  derivatives  based  upon  a  seimempirical  method 
due  to  Wood.  A  computer  program  formulation  of  Wood's  method  is 
presented. 
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SYMBOLS 

drag  coefficient,  D/QS 
normal-force  coefficient,  -Fz/QS 

normal-force  derivative,  3CN/3a 

pitching-moment  coefficient,  M^/QSd 
pitching-moment  derivative,  SC^/Sa 

damping-in-pitch  derivative,  3(^/3 (qd/2V)  +  3Cm/3 (ad/2V) 

rolling-moment  coefficient,  M,  /QSd 
damping-in-roll  derivative,  3C^/3 (pd/2V) 

yawing-moment  coefficient,  M  /QSd 

z 

Magnus  moment  coefficient,  3Cn/3 (pd/2V) 

Magnus  moment  derivative,  3 2Cn/3 (pd/2V) 3a 

side-force  coefficient,  F^/QS 
Magnus  force  coefficient,  30^/3 (pd/2V) 

Magnus  force  derivative,  3 2C^/d  (pd/2V) 3a 
reference  length 
drag  force 

force  along  the  Y  axis 
force  along  the  Z  axis 

moment  of  inertia  about  the  X  axis,  axial  moment 
of  inertia 

moment  of  inertia  about  the  Y  axis,  transverse  moment 
of  inertia 
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K 


x 

K 


y 


M 


x 


M 


nondimensional  axial  radius  of  gyration, 

nondimensional  transverse  radius  of  gyration, 

V^yy/md'7 

rolling  moment 
pitching  moment 


M2 

m 

P 

pd/2V ,  p 

Q 

S 

l/sg 

sd 

V 

V 

w 

x,y,z 

X 

y 

a 


*R 


R 


yawing  moment 
body  mass 

spin  rate  (rad/sec) 
reduced  spin  rate 
dynamic  pressure,  (l/2)  pV2 
reference  area 

reciprocal  of  the  gyroscopic  stability  factor 

damping  stability  factor 

airspeed 

vehicle  velocity  along  Y  axis 
vehicle  velocity  along  Z  axis 
conventional  aeroballistic  body  axes 
horizontal  coordinate  axis,  range 
vertical  coordinate  axis,  altitude 
angle  of  attack,  tan  ^  w/V 
total  angle  of  attack,  yfaJ  +  6 2 

angle  of  side  slip,  tan  ^  v/V 
yaw  angle  of  repose 

flight  path  angle 


3 


NOLTR  71-95 


TEST  FACILITIES 

All  aerodynamic  data  presented  in  this  report  were  measured 
in  the  Naval  Ordnance  Laboratory  Pressurized  Ballistics  Range 
No.  3  and  Supersonic  Tunnel  No.  1.  The  Pressurized  Ballistics 
Range  is  a  three-foot-diameter  170-foot-long  steel  tube  that 
can  be  evacuated  to  a  pressure  of  3mm  Hg.  There  are  27 
divergent  light  spark  shadowgraph  stations  along  the  range  tube 
spaced  alternately  five  and  eight  feet  apart.  At  these 
stations,  high-quality  shadowgraph  photographs  are  obtained 
in  both  the  vertical  and  horizontal  planes .  A  complete 
projectile  trajectory  is  determined  from  this  photographic 
information.  These  shadowgraphs  are  very  useful  for  flow 
phenomena  studies  also,  although  no  flow  visualization  was  used 
in  the  present  work. 

Model  launchers  available  for  use  in  this  facility  include 
a  variety  of  powder  guns,  both  rifled  and  smooth  bore,  as  well 
as  a  two-stage  light-gas  gun  with  a  smooth-bore  launch  tube 
of  20mm  diameter  and  a  smooth-bore  launch  tube  of  1.25-inch 
diameter.  In  these  tests  of  the  20mm  projectile  a  rifled 
powder  gun  was  used.  For  the  smooth  projectile  tests,  a  driver 
plug  was  used  to  engage  the  gun  rifling,  with  the  plug  in  turn 
keyed  to  the  projectile's  base. 

The  wind-tunnel  measurements  were  carried  out,  as  mentioned, 
in  the  NOL  Supersonic  Tunnel  No.  1.  This  is  a  blowdown  facility, 
having  a  16-  by  16-inch  test  section.  Stagnation  conditions 
are  essentially  atmospheric.  The  Pressurized  Ballistics  Range 
No.  3  and  Supersonic  Tunnel  No.  1,  as  well  as  other  NOL  Aero- 
and  Hydroballistics  Facilities,  are  described  and  illustrated 
in  Reference  (1) . 


CONFIGURATIONS 


All  aerodynamic  data  presented  in  this  report  are  centered 
around  the  Basic  20mm  General  Purpose  Projectile  or  configurations 
formed  from  simple  configurational  modifications  of  this  basic 
shape.  The  Basic  configuration  is  illustrated  in  Figure  1.  It 
will  be  noted  immediately  that  this  projectile  has  a  number  of 
surface  irregularities  due  to  nose  fuze,  rotating  band  and  crimping 
groove.  While  all  these  contour  undulations  are  necessary  on  the 
operational  projectile,  it  was  decided  to  construct  a  Smooth 
configuration  as  a  datum  to  ass-ess  the  effect  of  such  irregularities 
on  projectile  performance.  Thus,  the  Smooth  configuration  was 
constructed  from  the  Basic  shape  by  removing  the  irregularities. 

This  Smooth  configuration  is  illustrated  in  Figure  2.  Most  of  the 
data,  and  all  of  the  analysis,  in  this  report  are  concerned  with 
these  two  configurations. 

In  the  initial  phases  of  this  program  other  configurations  were 
studied.  Two  overlong  shapes  were  formed  by  adding  0.3-  and  0.6- 
caliber  cylindrical  afterbody  extensions  to  the  Basic  shape.  Both 
of  these  shapes  had  satisfactory  aerodynamic  performance,  as  will 
be  pointed  out  subsequently.  A  0.5-caliber  long,  seven-degree 
boattail  was  also  investigated. 

Some  of  the  relevant  physical  properties  of  the  Basic 
configuration  of  the  20mm  general  purpose  projectle  are  given  below: 


Length: 

Maximum  diameter: 

Center  of  gravity: 

Mass: 

Axial  moment  of  inertia: 
Transverse  moment  of  inertia: 

The  geometric  details  of  the  Basic 
Figure  3. 


4.9  calibers 

20  millimeters  or  0.786  inch 
2.21  calibers  from  base 
1.00578  x  10-2  slugs 
5.283  x  10"“  slug-ft2 
6.0133  x  10-3  slug-ft2 

configuration  are  given  in 


DATA  ACQUISITION 

In  the  Pressurized  Ballistics  Range,  spark  shadowgraphic 
photographs  are  taken  of  the  model  as  it  passes  each  station.  These 
photographs  are  etched  on  glass  to  provide  a  permanent  record.  Model 
attitude  and  down-range  position  are  read  on  an  optical  analyzer. 

Drag  and  aerodynamic  derivatives  are  then  reduced  from  these 
measurements  by  the  methods  outlined  in  Reference  (2) . 

In  the  wind-tunnel  tests  the  model  is  mounted  on  a  four-component 
strain-gage  balance.  Since  a  Magnus  test  requires  that  load 
measurements  be  made  with  the  model  spinning,  it  is  necessary  to 
provide  means  for  supporting  the  model  in  spin  to  keep  vibratory 
interaction  with  the  load-sensing  gages  to  a  minimum.  In  addition, 
there  must  be  some  source  of  torque  to  drive  the  model  in  spin. 


NOLTR  71-95 


Figure  4  illustrates  the  essentials  of  the  Magnus  balance. 

In  this  figure  a  representative  Magnus  model  is  shown  mounted  on  a 
Magnus  balance.  It  will  be  noted  that  the  model  makes  contact  with 
the  balance  via  fore  and  aft  ball  bearings.  Torque  is  generated 
by  means  of  a  two-stage  air  turbine  with  air  entering  through 
the  sting  and  exiting  from  the  model  base.  The  operational  procedure 
is  something  like  the  following:  Air  is  admitted  to  the  air  turbine 
at  about  400  psi.  Once  the  model  has  attained  a  spin  rate  of  about 
500  revolutions  per  second,  the  air  to  the  turbine  is  terminated. 

The  strain  gages  are  sampled  about  80  times  per  second  as  the  model 
undergoes  spin  decay.  A  magnetic  tachometer,  located  in  the  model, 
as  shown  in  Figure  4,  provides  a  nearly  continuous  record  of  spin 
rate.  If  gage  readings  are  converted  to  aerodynamic  coefficients 
by  means  of  the  balance  calibration  and  a  knowledge  of  the  flow 
properties,  it  is  easy  to  conceive  of  a  graphical  record  of  force 
or  moment  coefficient  versus  spin.  Actually,  such  a  record  is 
essential  to  wind-tunnel  Magnus  data  reduction. 

The  view  indicated  in  Figure  4  is  such  that  the  plane  of  the 
paper  is  the  yaw  plane.  The  component  indicated  as  the  Magnus 
flexure  lies  at  the  heart  of  the  Magnus  balance.  A  set  of  these 
Magnus  flexures  is  illustrated  in  Figure  5.  It  will  be  noted  that 
the  flexure  is  essentially  a  miniature  beam  mounted  such  that  it 
will  be  loaded  eccentrically  by  the  main  beam  or  balance  as  the 
balance  is  subject  to  bending  by  the  aerodynamic  loads  on  the  model. 
For  example,  as  the  balance  is  subjected  to  yaw  loads  the  flexure 
is  eccentrically  end  loaded.  This  induced  secondary  bending  results 
in  the  flexure  acting  as  a  mechanical  amplifier.  This  secondary 
bending  is  illustrated,  although  with  considerable  exaggeration, 
in  Figure  6 ,  Because  the  Magnus  flexures  are  used  only  in  the 
yaw  plane,  the  sensitivity  of  the  balance  in  yaw  is  about  nine 
times  that  of  the  balance  in  pitch. 

DATA  REDUCTION 

No  discussion  will  be  given  herein  of  the  methods  used  in 
reducing  ballistics  range  measurements.  An  outline  of  the 
mathematical  foundations  is  given  in  Reference  (2) .  Several 
formulations  of  ballistics  range  data-reduction  methods  have  been 
made  for  application  to  the  technique  of  wind-tunnel  free-f light. 

One  such  report  is  Reference  (3) . 

In  the  wind  tunnel,  force  and  moment  data  were  obtained  using 
the  Magnus  balance  in  a  rotating  sector  arm.  The  model-balance 
combination  is  brought  to  a  desired  angle  of  attack,  the  model  then 
spun  to  the  desired  spin  rate  and  the  drive  air  terminated  after 
the  flow  conditions  are  established.  After  tunnel  flow  has  been 
established,  the  tachometer  and  the  strain  gages  are  sampled  at 
about  80  times  per  second.  These  signals  are  digitized  and 
recorded  on  magnetic  tape.  This  record,  together  with  the 
balance  calibration  and  flow  characteristics,  produces  plots  of  the 
force  and  moment  coefficients  as  a  function  of  reduced  spin  rate  at 
a  specific  value  of  Mach  number  and  angle  of  attack.  The  variation 
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of  normal  force  and  pitching  moment  with  spin  is  usually  small  and 
considered  negligible  in  this  presentation.  However,  the  side-force 
and  yawing-moment  coefficient  indicate  a  nearly  linear  variation 
with  spin  rate.  A  representative  plot  of  side-force  and  yawing- 
moment  coefficients  versus  reduced  spin  rate  is  given  in  Figure  7. 

It  will  be  noted  that  even  at  this  rather  large  angle  of  attack 
these  force  and  moment  coefficients  are  linear  with  reduced  spin 
rate.  In.  order  to  discuss  the  Magnus  data-reduction  methods  as 
used  in  wind  tunnels,  it  is  necessary  to  consider  the  analytic 
formulation  of  the  Magnus  effect. 

The  Magnus  force,  F,  will  be  defined  as  a  force  depending  upon 
body  spin  rate  and  angle  of  attack,  and  acting  normal  to  the  plane 
established  by  the  spin  vector,  p,  and  £he  free-stream  velocity 
vector,  V^.  Mathematically  this  force,  F,  and  its  corresponding 

moment,  can  be  expressed  as 

f -*(!»■*  \L  <l> 


M  -  k  [f  *  &] 


where  k  is  a  scalar  constant  for  a  given  set  of  flow  conditions; 
r  is  the  vector  distance  from  the  center  of  gravity  to  the  Magnus 
center  of  pressure,  along  the  body's  axis  of  symmetry;  and,  p  i,s 
the  spin-rate  vector  defined  along  the  axis  of  symmetry. 

All  forces  and  moments  are  referred  to  the  conventional  aero- 

ballistic  body  axis  system,  an  axis  system  which  is  fixed  to  the  body 

and  shows  all  its  rotational  and  translational  motion,  except  spin. 

In  the  axis  system,  the  x  axis  is  forward  along  the  axis  of  symmetry; 

the  y  axis  is  to  the  right  when  the  store  is  viewed  along  the 

positive  x  axis;  the  z  axis  completes  a  right-handed  triad.  The 

origin  of  this  axis  system  is  at  the  moment  reference  center  taken, 

in  this  case,  to  be  the  center  of  gravity^.  Unit  vectors  along  the 

x,  y,  z  axes  will  be  defined  as  1  ,  1  ,  1  .  The  wind-tunnel 
J  x  y  z 

constraints  are  such  that  the  x,  z  plane  is  vertical;  and,  that 
this  plane  contains  the  flow  velocity  vector. 

Since  p  =  p  1  ,  V  =  (V  cosa)l  +  (V  sina)t  ,  Equation 

may  be  rewritten  as 


p  —  -  K  V  &  Fyi-y 
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where  the  side  force,  F  ,  is  equal  to  -kpV  sina.  This  relationship 

y 

demonstrates  that  the  Magnus  force  is  an  odd  function  of  the  angle 
of  attack  with  spin  rate.  The  Magnus  moment  equation  also  may  be 
rewritten  using  the  components  of  the  spin-rate  and  free-stream 
velocity  vectors.  That  is, 

M- -K'LS’rSin*]!,'  (,1> 

where  the  yawing  moment,  M  ,  is  equal  to  -kpr  sina.  It  can  be 
seen  from  Equation  (4)  that  the  Magnus  moment  is  an  odd  function  in 
center-of-pressure  location,  spin  rate  and  angle  of  attack.  For 
example,  the  Magnus  center  of  pressure  is  usually  behind  the  center 
of  gravity,  giving  a  positive  Magnus  moment  (nose  to  the  right  as 
viewed  from  the  rear) . 

The  yaw-force  and  vaw-moment  coefficients  are  defined  as: 


<?r=  Fy/QS  Cn=  Mt/QSd  (-5) 

The  above  coefficients  depend  upon  the  body  pressure  distribution, 
which,  in  turn,  depends  upon  the  compressibility,  viscosity  and 
unsteadiness  of  the  flow  field.  To  indicate  the  degree  of  simulation 
of-  these  effects,  it  is.  necessary  to  present  coefficients  as  functions 
of  the  appropriate  flow  similarity  parameters. 

Since  the  free-stream  velocity  is  in  the  vicinity  -of  the 
speed  of  sound,  it  is  necessary  to  regard  the  medium  as  compressible. 
Simulation  of  compressibility  effects  is  assured  by  testing  at 
identical  free-f light  Mach  numbers.  Also,  since  the  Magnus  effect 
on  the  forebody  originates  entirely  in  the  boundary  layer,  it  is 
necessary  to  test  at  the  anticipated  Reynolds  numbers  to  simulate 
viscous  effects.  Finally,  since  each  surface  element  on  a  steadily 
spinning  body  experiences  a  cyclically  changing  flow  field,  the  test 
must  be  made  at  a  parameter  which  matches  flow  unsteadiness.  In 
Magnus  tests  this  flow  unsteadiness  parameter  is  designated  as  the 
reduced  frequency,  p.  Testing  at  identical  reduced  frequencies 
assures  a  matching  of  the  flow  angularity  at  similarly  located 
surface  elements  on  geometrically  similar  bodies.  Thus,  it  rill 
be  postulated  in  a  Magnus  test  that  the  coefficients  must  be 
expressed  as  functions  of  Mach  number,  Reynolds  number  and  reduced 
frequency,  as  well  as  body  angular  attitude. 


If  the  Magnus  force  is  assumed  to  be  an  analytic  function  of 
angle  of  attack  and  reduced  frequency,  the  yaw-force  coefficient  can 
be  expanded  in  a  truncated  Taylor  series  in  a  and  p,  as: 


-v  ^ +  + tf  f 
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where  all  derivatives  are  evaluated  at  a  and  p  equal  to  zero.  Since, 


G  (<y)  ’  V'M*0 


(7) 


it  follows  that  all  but  the  cross  derivatives  vanish.  Thus,  Equation 
(6)  becomes,  as  a  first  approximation, 


(8) 


A  similar  relationship  for  the  moment  coefficient,  C  ,  would  be 


^r*f  ■ 

The  terms  on  the  right  in  Equations  (8)  and  (9)  are  the 
familiar  Magnus  force  and  moment  derivatives,  respectively,  for  linear 
aerodynamics.  It  should  be  noted  that  Equations  (8)  and  (9)  are 
compatible  with  Equations  (3)  and  (4)  for  small  angles  of  attack. 


Before  coming  to  grips  with  an  analytic  description  of  the 
problem's  nonliriearities,  it  is  important  to  recall  the  methods 
of  data  acquisition.  Side-force  and  yawing-moment  measurements  were 
made  while  the  model  was  undergoing  spin  decay  at  a  fixed  angle 
of  attack  and  immersed  in  a  flow  of  constant  Mach  number.  After 
tachometer  and  strain-gage  signals  are  recorded,  the  angle  of  attack 
is  changed.  The  model  is  spun  again  to  some  predetermined  upper 
limit,  the  drive  air  is  terminated  and  the  procedure  repeated. 

After  the  entire  angle-of-attack  range  has  been  spanned,  the  Mach 
number  is  changed  and  the  procedure  repeated.  Since  it  is  noted 
that  for  configurations  without  fins,  the  Magnus  force  and  moment 
coefficients  are  linear  with  spin  rate,  Equation  (6)  may  be 
rewritten  as 


^CyU)/3jp=  (10) 

The  variation  of  the  Magnus  effect  with  angle  of  attack  is  represented 
by  the  derivative,  3Cy(a)/9p.  If  the  Magnus  effect  varies  linearly 

with  angle  of  attack,  then  Equation  (10)  may  be  replaced  by  Equation 
(8) .  Experiments  show  that,  generally,  the  Magnus  effect  varies  in 
a  nonlinear  fashion  for  angles  of  attack  greater  than  two  or  three 
degrees.  Thus,  the  wind-tunnel  procedure  is  to  obtain  OC^/Op  for 

several  angles  of  attack  at  a  fixed  Mach  number.  Since  the 
ballistics  range  data-reduction  methods  are  based  upon  linear 
aerodynamics,  the  Magnus  effect  can  only  be  represented  by  the  Magnus 
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derivatives  of  Equations  (8)  and  (9).  Such  a  limitation  is  usually 
acceptable  as  ballistics  range  tests  are  usually  limited  to 
a  few  degrees  angle  of  attack.  Wind-tunnel  measurements  can  be 
compared  with  ballistics  range  measurements  by  obtaining  the  slope 
through  the  origin  of  the  SC^/3p  versus  angle-of-attack  graphs. 

This  slope,  3/3a(3C  /3p)  or  C  is  the  Magnus  force  derivative. 

^  ^pa 

AERODYNAMIC  MEASUREMENTS 

As  mentioned  several  times  earlier,  the  aerodynamic  data 
contained  in  this  report  were  obtained  in  the  NOL  Pressurized 
Ballistics  Range  and  Supersonic  Tunnel  No.  1.  The  flow  environment 
experienced  by  the  range  and  wind-tunnel  models  might  be  summarized 
in  a  plot  of  Reynolds  number  versus  Mach  number.  Such  a  plot  is 
presented  in  Figure  8  with  the  Reynolds  number  being  based  upon 
body  diameter. 

In  the  ballistics  range,  two  pressurization  levels  were  used — 
one  atmosphere  and  one-tenth  of  an  atmosphere.  These  levels  are 
also  indicated  in*  Figure  8.  The  wind  tunnel,  being  a  blowdown 
facility,  has  a  total  pressure  head  of  one  atmosphere.  At  the  lower 
Mach  numbers  (below  Mach  1.5)  it  will  be  noted  that  a  higher  Reynolds 
number  is  achieved  in  the  wind  tunnel.  This  higher  Reynolds  number 
occurs  because  the  wind-tunnel  model  is  nearly  three  times  the 
linear  size  of  the  full-scale  ballistics  range  models  (2.0  inches  as 
compared  to  0.786  inch  for  the  ballistics  range  model). 

Figures  9  through  17  present  the  wind-tunnel  measurements  of 
the  normal-force  and  pitching-moment  coefficients  for  the  Basic 
configuration;  Figures  18  through  24  present  the  same  information 
for  the  Smooth  configuration.  It  will  be  noted  that  the  pitching- 
moment  coefficient  is  linear  through  the  origin  and  remains  linear 
up  to  about  10  degrees  angle  of  attack  fcr  subsonic  measurements 
(see,  for  example,  Fig.  11) .  Supersonically,  the  pitching  moment 
begins  to  give  evidence  of  nonlinearities  at  about  eight  degrees 
(see  Fig.  15) .  This  same  sort  of  observation  can  be  made  for  the 
Smooth  body,  indicating  that  nose  geometry  rar.her  than  surface 
irregularities  is  the  main  contributor. 

It  was  mentioned  earlier  that  a  variety  of  configurations  were 
briefly  examined  in  the  ballistics  range  at  the  inception  of  this 
program.  In  addition  to  the  Smooth  conf inu'-ation,  these  configurations 
consisted  of  a  0.3-caliber  cylindrical  e  ;ion,  a  0.6-caliber 
cylindrical  extension  and  a  seven-degre  . ^-caliber-long 
boattail.  The  normal-force  and  pitchir  w  A  it  derivatives  for 

these  configurations  are  compared  with  x  .  ■  of  the  Basic  body  in 
Figures  25  and  26.  It  will  be  noted,  for  example,  that  the  0.3- 
caliber  extension  gives  a  pitching  moment  about  equal  to  that  of  the 
Basic  configuration;  the  0.6-caliber  extension  seems  to  make  an 
appreciable  increase  in  the  pitching  moment.  The  boattail  configuration 
has  a  significantly  lower  pitching  moment.  It  should  be  emphasized 
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that  the  Smooth  body  data,  given  in  Figures  25  and  26,  are 
questionable  and  should  not  be  used  for  comparative  purposes.  Due 
to  a  defective  plug  design  there  was  gas  blow-by  which  seemed  to 
distort  the  model  base.  This  design  was  corrected  in  later  tests. 

Since  the  most  important  configuration  is  these  tests  is  the 
Basic  body,  it  was  desired  to  have  a  fairly  detailed  data  coverage. 
Figures  27  and  28  present  comparisons  of  ballistics  range  and  wind- 
tunnel  measurements  of  the  normal-force  and  pitching -moment  derivative's 
of  the  Basic  body.  Figure  29  gives  the  center-of-pressure  location 
from  the  body  vertex,  calculated  from  the  data  of  Figures  27  and  28. 

The  correlation  between  ballistics  range  and  wind-tunnel  data  is 
quite  satisfactory.  It  should  be  recalled  from  Figure  8,  that 
the  Reynolds  numbers  in  both  facilities  are  fairly  close  up  to  a 
Mach  number  of  1.5.  Figure  29  shows  that  the  center  of  pressure 
moves  rearward  with  increasing  Mach  number.  Figures  30  and  31 
present  the  normal-force  and  pitching -moment  derivatives,  respectively, 
for  the  Smooth  configuration.  The  data  for  the  Smooth  configuration 
are  much  more  sparse  than  the  Basic  configuration,  with  only  one 
data  point  available  from  the  ballistics  range.  Figure  32  presents 
the  center-of-pressure  location  in  calibers  aft  of  the  vertex.  In 
comparing  Figure  29  to  32  it  appears  that  even  within  the  scatter 
of  the  data,  the  center  of  pressure  for  the  Smooth  configuration 
is  slightly  ahead  of  that  for  the  Basic  configuration. 

One  of  the  most  important  single,  measurement  made  in  the 
ballistics  range  is  drag.  Drag  measurements  for  the  Basic  and 
Smooth  configurations  are  presented  in  Figure  33  in  the  form  of 
the  zero-lift  drag  coefficient.  It  is  quite  obvious  from  this 
figure  that  the  Smooth  configuration  enjoys  a  significant  drag 
advantage  over  the  Basic  configuration. 

In  Figures  34  and  35  the  damping-in-roll  derivative,  , 

P 

and  the  damping-in-pitch  derivative,  C  +  C  ,  are  presented 

q  a 

respectively,  as  functions  of  Mach  number.  The  damping-in-roll 
derivative  is  reduced  according  to  the  methods  of  Reference  (4) . 

The  roll  angle  is  determined  by  pins  located  in  the  base  of  the  model. 
The  damping-in-pitch  derivative,  unfortunately,  has  a  great  deal 
of  scatter  in  the  measurements.  Since  this  derivative  is  essential 
in  specifying  the  dynamic  stability  criterion,  S^,  there  will  be 

some  uncertainty  concerning  regions  of  dynamic  stability. 

Also  vital  in  the  determination  of  the  dynamic  stability  of 
a  configuration  is  the  Magnus  effect.  The  Magnus  effect  had  been 
discussed  earlier  in  the  Data-Reduction  Section.  Magnus  measure¬ 
ments  were  carried  out  in  both  the  wind  tunnel  and  ballistics  range. 

In  the  wind  tunnel,  both  the  side  force  and  yawing  moment  are 
measured  as  functions  of  spin  rate  and  angle  of  attack.  As  pointed 
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out  in  Equation  (10) ,  the  force  and  moment  are  considered  linear  with 
spin  rate  and  are  expressible  as  force  and  moment  derivatives 
with  spin  rate.  These  derivatives,.  SC^/Sp  and  3Cn/3p,  have  been 

designated  as  the  Magnus  force  and  moment  coefficients,  C  and  C  , 

yp  p 

respectively.  Figures  36  through  41  present  the  Magnus  force  and 
moment  coefficients  for  the  Basic  configuration  and  Figures  .42  through 
47  for  the  Smooth  configuration.  It  will  be  noted  that  while  the 
slope  of  the  Magnus  force  coefficient,  the  Magnus  force  derivative, 
remains  negative  over  the  entire  Mach,  number  range,  the  Magnus  moment 
derivative  changes  sign  in  going  from  subsonic  to,  supersonic  Mach 
numbers  (compare  Figure  37  to  38) .  The  above  observation  indicates 
that  the  Magnus  center  of  pressure  is  slightly  ahead  of  the  center 
of  gravity  subsohically  and  slightly  aft  of  the  center  of  gravity 
supersonically .  it  will  be  pointed  out  later  that  a  Magnus  center 
of  pressure  ahead  of  the  center  of  gravity  always  results  in  dynamic 
instability.  Thus,  it  would  be  expected  that  the  Basic  configuration 
is  dynamically  unstable  subsonically . 

Interestingly  enough,  the  Smooth  configuration  appears  to  be 
dynamically  stable  over  the  whole  flight  regime  as  the  slope  of  the 
Magnus  moment  coefficient,  Magnus  moment  derivative,  is  negative  for 
subsonic  and  supersonic  Mach  numbers  (compare  Figure  43  to  44) .  The 
dynamic  instability  of  the  Basic  configuration  is  not  serious  as  it 
would  be  expected  that  this  projectile  would  only  be  effective  in 
the  moderately  high  Mach  number  range  (Mach  2  to  3.5) . 

In  examining  the  Magnus  data  of  Figures.  36  through  47  it  will 

be  noted  that  the  variation  of  C  and  C  with  angle  of  attack  is 

y~  n 

P  P 

nonlinear  for  angles  much  beyond  two  degrees..  This  is  particularly 
true  for  the  Basic  configuration.  For  small  angles  of  attack  there 
appears  some  justification  in  linearizing  the  Magnus  effect  by 
representing  the  Magnus  force  and  moment  by  the  Magnus  derivative 
(see  Eq.  (10)).  This  linearization  has  been  done  for  the  Basic 
configuration  and  the  resulting  Magnus  force  and  moment  derivatives 
plotted  as  functions  of  Mach  number  in  Figure  48.  These  data  appear 
quite  ragged  because  the  Magnus  derivative  is  formed  by  twice 
differentiating  experimental  data,  once  with  respect  to  spin  rate 
and  again  with  respect  to  angle  of  attack. 

In  the  ballistics  range  an  entirely  different  procedure  is 
used.  Since  the  observed  motion  (as  obtained  from  the  photographic 
plates)  is  assumed  to  be  the  "solution"  of  a  set  of  linear 
differential  equations,  the  Magnus  effect  can  only  be  expressed  in 
terms  of  the  Magnus  derivative.  Also,  as  measurements  were  made  at 
only  one  center-of -gravity  location,  only  the  Magnus  moment 
derivative  is  available. 

Figure  49  presents  the  Magnus  moment  derivative  as  a  function 
of  Mach  number  comparing  the  wind-tunnel  and  ballistics  range 
measurements.  It  will  be  observed  that  the  Magnus  moment  changes 
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sign  somewhere  around  a  Mach  number  of  1.1.  Thus,  one  would  expect 
that  the  projectile  would  become  dynamically  unstable  below  this 
Mach  number. 

I  ANALYTIC  DETERMINATION  OF  THE  STATIC  COEFFICIENTS 

;  There  are  a  variety  of  methods  for  predicting  the  normal-force 

|  and  pitching-moment  derivatives  as  a  function  of  Mach  number. 

■,  '  Rrobably  the  simplest  is  that  due  to  Hitchcock  (Ref.  (5)). 

\  Hitchcock's  method  is  entirely  empirical  and  has  the  disadvantage 

\  of  ignoring  compressibility  effects.  The  normal-force  derivative, 

l  '  CN  '  an<^  Center  pressure,  h,  are  given  as  follows: 

?  a 

i 

4 

!  Cw=  |:[.(!,53+.oZ23l\~.fci3^B-.oo^C+.3.635‘D+-^4730/E')l  (11) 


(12) 


where 

A  is  the  angle  of  the  boattail  in  degrees 
B  is  the  length  of  the  boattail  in  calibers 

C  is  the  length  of  the  cylindrical  part  in  calibers 

D  is  the  length  of  the  ogival  head  in  calibers 

E  is  the  radius  of  the  ogival  arc  in  calibers 

If  reference  Is  made  to  Figure  27  it  is  possible  to  compare  the 
prediction  of  CN  ,  using  Hitchcock's  formula  of  Equation  (11),  with 

a 

wind-tunnel  and  ballistics  range  measurements. 


With 

A  =  0 
B  =  0 
C  =  2.41 
D  =  2.49 
E  =  7.40 


Equation  (11)  gives 


3.54 
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This  value,  which  is  invarient  with  Mach  number,  is  higher  than 

any  measurement  made  in  either  the  wind  tunnel  or  ballistics  range 

at  any  Mach  number.  According  to  Equation  (12)  the  center  of 

pressure  is  2.17  calibers  aft  of  the  vertex.  The  Hitchcock 

relationship  predicts  a  center  of  pressure  which  is  too  far  aft. 

The  pitching-moment  derivative,  C  ,  determined  from  Equations  (11) 

m 

a 

and  (12),  is  1.85.  Comparison  with  Figure  28  shows  that  this  value 
is  lower  then  any  of  the  measurements,  regardless  of  Mach  number. 
Supposedly,  the  empirical  constants  of  Equations  (11)  and  (12)  could 
be  adjusted  to  bring  the  predictions  in  better  agreement  with  the 
measurements.  However,  the  primary  shortcoming  of  this  method  is, 
as  pointed  out  earlier,  its  failure  to  account  for  compressibility 
effects. 

A  second  method  due  to  Wood  (Ref.  (6))  is  available  which 
does  take  into  account  compressibility  effects.  Wood  defines  the 
normal-force  derivative  and  the  pitching-moment  derivative  about 
the  base,  respectively,  as: 

[(ZSt/S>.5]f  (13) 


where  =  (S  +  S^)/2  and  where  S  and  are  the  reference  area  and 

base  area,  respectively.  The  center-of-pressure  location,  in 
calibers,  forward  from  the  base  is  given  as, 


c 


C 


b 


■H  =  { 

The  pitching-moment  derivative  about  the  center  of  gravity  is 
easily  written  from  Equations  (14)  and  (15)  as, 


CM  -  GjC^ 


(16) 


where  G  is  the  location  of  the  center  of  gravity,  in  calibers, 
forward  of  the  base.  The  functions,  f^  and  f£,  are  empirical 

functions  of  the  fineness  ratio  and  the  Mach  number.  These 
functions  are  equal  to  unity  for  Mach  numbers  between  0  and  unity 
and  are.  available,  in  graphical  form,  for  Mach  numbers  in  excess 
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of  unity.  Graphs  of  f^  and  f£  are  available  in  both  References 
(6)  and  (7)  . 

Table  1  presents  a  listing  of  a  BASIC  program  to  calculate 
the  normal-force  and  pitching-moment  derivatives  and  center-of- 
pressure  location  aft  of  the  vertex.  In  order  to  use  the  f^  and 

f 2  functions  in  the  program,  a  seventh-order  polynomial  was 

fitted  to  these  graphs.  These  polynomials,  together  with  Equations 

(13)  through  (16) ,  were  formulated  in  the  BASIC  language  to 

provide  a  fairly  good  semiemphirical  estimate  of  the  normal-force 

and  pitching-moment  derivatives  and  the  center  of  pressure. 

Estimates  of  CM  ,  C  and  center-of-pressure  location  from  Wood's 
N  m  c 

a  a 

method  are  compared  with  wind-tunnel  and  ballistics  range  measurements 
in  Figures  27,  28  and  29.  The  agreement  may  be  considered  fair 
and  probably  limited  to  preliminary  design  estimates. 

BASIC  COMPUTER  PROGRAM  FOR  THE  RAPID  DETERMINATION 
OF  PROJECTILE  STABILITY 

In  this  section,  a  computer  program  developed  for  the  rapid 
evaluation  of  the  stability  of  an  unguided  weapon  is  described. 

This  program  assuir.es  that  the  forces  and  moments  acting  on  the 
body  arise  from  linear  aerodynamic  loads  and  a  constant  gravitational 
field.  More  elaborate  simulation  of  projectile  characteristics 
is  certainly  possible  in  six-degree-of-freedom  programs,  such  as 
Reference  (8) .  While  the  increased  precision  of  such  programs  is 
obvious,  the  lengthy  table  preparation  and  long  computation  time  often 
discourages  the  use  of  such  programs.  For  example,  it  is  not 
uncommon  to  have  a  computation  time  2000  times  longer  than  real 
time.  In  this  program  only  two  second-order  differential  equations 
are  solved.  Other  relationships,  which  might  be  written  in  the 
form  of  differential  equations,  are  approximately  integrated  and 
it  is  these  analytic  approximations  that  are  programmed.  The 
assumption  of  linear  aerodynamics  avoids  lengthy  tables.  This 
latter  assumption  may  not  be  severely  restricting,  as  a  designer's 
first  performance  studies  are  usually  limited  to  estimates  of 
the  stability  derivatives.  Even  in  considering  operational 
vehicles,  satisfactory  performance  usually  restricts  the  projectile 
to  small  angular  excursions. 

This  program  was  written,  primarily,  for  spin-stabilized 
projectiles  although  it  is  entirely  adequate,  within  the  restrictions 
of  the  force  and  moment  representations,  to  spinning  fin-stabilized 
weapons  as  well.  Regardless  of  the  stabilization  characteristics 
of  the  weapon,  the  program  operates  in  the  following  manner: 

A  particle  trajectory  is  calculated  to  provide  the  dynamic  pressure, 
velocity  and  Mach  number.  These  quantities  are  then  used  to 
calculate  the  reduced  frequency,  the  yaw  of  repose,  the  damping 
stability  parameter  and  the  gyroscopic  stability  parameter. 
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Various  print  statements  provide  these  quantities  as  a  tabular 
function  of  altitude  and  time. 

In  carrying  out  the  computation  there  are  three  distinct 
sets  of  equations.  First,  two  second-order  differential  equations, 
nonlinearily  coupled  through  aerodynamic  drag,  are  used  to  provide 
a  particle  trajectory  program.  This  particle  trajectory  "locates" 
the  weapon  as  it  moves  through  a  uniform  gravitational  field 
subject  to  aerodynamic  drag.  Since  a  solution  in  closed-form 
is  out  of  the  question,  it  was  necessary  to  solve  these  equations 
using  numerical  integration. 

The  second  set  of  equations  provides  reduced  frequency  as  a 
function  of  arc  length.  If  the  Mach  number  is  assumed  to  be 
constant,  this  equation  may  be  integrated  to  give  reduced  frequency 
as  a  function  of  arc  length.  If  changes  in  Mach  number  affect 
the  value  of  the  drag  and  roll-damping  coefficients,  it  is  necessary 
to  carry  out  the  integration  piecewise  over  intervals  of  constant 
Mach  number.  Such  a  procedure  avoids  a  lengthy  numerical  integration 
procedure  and  keeps  the  program  formulation  within  a  reasonable 
length. 

Finally,  a  third  set  of  equations  is  needed  to  describe  the 
angular  motion  of  the  body.  Two  coupled  second-order  differential 
equations  are  integrated  to  express  the  ~aw  and  pitch  motion  of 
the  vehicle  about  its  velocity  vector.  _ne  integrals  of  these 
equations  are  used  to  formulate  the  gyroscopic  and  damping 
stability  parameters.  As  in  the  case  of  the  spin-rate  equation, 
this  integral  of  the  angular  motion  is  valid  only  over  intervals 
of  constant  Mach  number.  The  result  is  that  the  above-mentioned 
stability  criteria  will  vary  with  Mach  number. 

The  sketch  given  below  illustrates  the  variables  used  in 
the  trajectory  equations.  It  will  be  assumed  that  all  motion 
takes  place  in  a  vertical  plane,  with  the  "y"  axis  vertical 
(defined  positive  upwards)  and  the  "x"  axis  horizontal  (position 
forward  in  the  directions  of  flight) . 

V 
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If  the  gravitational  and  drag  forces  are  resolved  into  the 
vertical  and  horizontal  components  (according  to  the  above 
sketch)  the  equation  of  vertical  and  horizontal  motion  may  be 
written  as, 


jVx  _  _  D  Qs&  _  _ 

in  mV 

Vx  =  d*  / dt 

dVy„  _  _ g  „ EL  -ft 

dt  mo  m  v  o 

Vy  =  Jy/dr 


(17a 

(17b 

(17c 

(17d 


Where  Equations  (17)  have  been  rewritten  in  alternate  form  by 

replacing  cos0  by  V  /V  and  sin6  by  V  /V.  If  the  following 

x  y 

substitutions  are  made: 


D  =  Al 

X 

=  X 

m  =  W 

dx 

=  XI 

dt 

g  =  32.174 

Vx 

=  Q 

V  =  V 


=  Q1 


y  =  Y 

&■.  VI 
dt 


Vy  P 

dV 

x  =  PI 
dt 


Equations  (17)  may  now  be  expressed  in  terms  of  computer  variables 
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All  computer  variables  are  defined  in  terms  of  problem  variables  in 
Table  2. 


The  second  equation  needed  in  the  stability  analysis  is  that 
relating  spin  rate  to  some  independent  variable,  such  as  time 
or  range.  The  most  straightforward  approach  is  to  first  state 
the  following  easily  established  identity. 


d  (  P\  iif .  t  dv 
dtvvj  vdt  v^dt 


(19) 


where  p  is  the  spin  rate,  t  is  the  time  and  V  is  the  free-stream 
airspeed.  From  the  roll-damping  equation,  or  the  equation  relating 
deceleration  to  the  roll-damping  moment,  we  have, 


db  _  1 

>QSd 

(  M  ^ 

i  _  ( cyv  s  A  b 

(20) 

ft 

V  i 

V  '  4m  r 

where  I  is  the  moment  of  inertia  about  the  longitudinal  axis,  C„ 
xx 


the  damping-in-roll  derivative,  m  the  mass  and  K  the  nondimensional 

X 

radius  of  gyration  defined  as  I  /m.  Next,  the  drag  equation  can  be 

XX 

written  as 


Inserting  Equations  (20)  and  (21)  into  identity,  (19)  gives 


The  above  equation  may  be  rewritten  changing  the  independent  variable 
from  time,  t,  to  downrange  distance,  x,  through  the  following 
operator 
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Equation  (22)  then  becomes. 


Equation  (23)  can  be  rearranged  to  give, 

<jQ>/v)  f  1 

t>/v  VCosG- 

where  cos9  =  V  /V.  An  exact  integral  of  Equation  (24)  in  closed- 

form  is  impossible  because  cos6  is  a  function  of  x,  which  can  only 
be  known  by  integrating  the  particle  trajectory.  However,  if 
the  integration  is  carried  out  over  a  short  enough  time  interval 
that  cos0  may  be  replaced  by  a  constant  value,  say  its  value  at 
the  beginning  of  an  interval  "i,"  the  integral  of  Equation  (24)  may 
be  written  as, 


If  the  number  5  is  attached  to  a  symbol  to  mean  the  value  at  the 
beginning  of  an  interval,  it  is  possible  to  rewrite  Equation  (25) 
in  computer  variables. 


(26) 

£5“  HXFCC V5YQS>C(CR*S)/a*W^  *  (Kk*lD *  .S  *  E  *  ft-  *$) 


P7=  CPSfcVVVS’ 

Now  P5,  the  value  at  the  beginning  of  the  next  interval,  is 

P5  =  P7*E5.  Again,  reference  is  made  to  Table  2  for  the  definition 

of  the  computer  variables  in  terms  of  the  problem  variables. 

Equations  (17)  describe  the  position  of  the  projectile  as 
a  function  of  time  as  it  moves  through  a  uniform  gravity  field. 

It  is  obvious  that  no  interaction  is  assumed  between  the  angular 
motion  of  the  body  and  its  trajectory.  This  is  obviously  a 
simplifying  assumption,  but  is  justified  on  the  grounds  that  the 
angular  motion  of  a  stable  projectile  should  be  small  and,  therefore 
interaction  terms,  such  as  lift  force,  Magnus  force  and  angular 
dependent  drag,  can  be  neglected.  The  question  arises  as  to  how 
one  can  be  sure  that  the  angular  motion  will  remain  small  under 
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response  to  an  initial  disturbance.  In  order  to  answer  this 
question,  two  stability  criteria  will  be  developed  to  provide 
the  necessary  and  sufficient  conditions  that  the  projectile's 
angular  motion  will  not  grow  in  response  to  a  disturbance.  These 
criteria  will  be  formulated  in  terms  of  the  body's  mass  distribution, 
spin  rate  and  Magnus,  damping  and  static  aerodynamic  characteristics. 

Murphy,  in  Reference  (2)  Equations  4.24  and  4.25,  introduces 
the  complex  variable 


c  y  4-iW 

c  v 

which,  for  small  angles  of  attack  and  sideslip,  becomes. 


(27a) 


^  +  i  &C  (27b) 

The  above  formalism  allows  the  differential  equations  of  lateral 
motion  to  be  expressed  in  terms  of  a  single  second-order  differential 
equation,  rather  than  in  two  second-order  equations  in  the  real 
variables  v/V  and  w/V  or  0  and  a.  As  indicated  by  Murphy  in  Equation 
6.12  of  Reference  (2),  this  equation  of  motion  becomes 


+  (H+i  p)£-  cm-h  rr)£  -  G 


(28) 


where 


(29a) 

(29b) 

(29c) 


(29c) 


Equations  (29a)  and  (29c)  differ  from  those  given  in  Reference  (2) 
(see  Eq.  6.12),  in  that  the  coefficients  C  +  C  and  Cm  are 

q  a  pa 

multiplied  by  1/2.  This  factor  of  1/2  is  necessary  as  NOL  dynamic 
moments  are  nondimensionalized  with  respect  to  pd/2V  whereas  in 
Reference  (2),  pd/V  is  used.  This  factor  of  1/2  occurred  earlier  in 
Equation  (23) . 
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Since  Equation  (28)  is  linear,  its  solution  is  given  as, 


Since  the  requirement  for  stability  is  that 


UW)  k  6 

S->oO 


it  is  necessary  and  sufficient  that  A^  and  X2  are  both  negative 
The  exponentials  of  Equation  (30)  may  be  written  as. 


\  v + i  4;,  -  i  Li  P-  H  i  J  4  M+H1-  v +  z  I  p  ( 2T-  H ) 

If  the  Magnus-lift  coupling  term,  T,  and  the  damping  term,  H, 
are  neglected,  one  has 

-  |[i  £  ill-  Vsjj] 


where  1/Sg  =  4M/P2  and  is  designated  as  the  reciprocal  of  the 

gyroscopic  stability  parameter.  Thus,  in  a  spin-stabilized  vehicle 
it  is  necessary  that  1/Sg  be  less  than  unity  for  the  vehicle  to 

be  stable.  A  fin-stabilized  vehicle  will  always  be  gyroscopically 

stable  as  M  is  negative.  If  1/S  is  greater  than  unity,  A  would 

9 

equal  1/yj 1/S  -  1  and,  being  positive,  would  result  in  divergence 

of  the  projectile  from  its  initial  conditions. 

If  1/S  is  less  than  unity  and  A,  and  A-  are  both  zero, 
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Equation  (33)  followed  from  Equation  (31)  after  T  and  H  were 
set  equal  to  zero.  The  exact  values  of  the  frequencies  ^  and 

will  vary  with  the  relative  magnitudes  of  T,  H,  M  and  P;  however. 
Equation  (33)  is  usually  an  excellent  approximation  for  the  modal 
frequencies  4*^  and  .  The  criteria  for  such  a  approximation  is 
that:  1  1 


Murphy  (Ref. 
values  of 


(2))  then  goes  on  to  show  that,  regardless  of  the 
t  e  A^  and  A2,  one  may  write, 


t,a. 


+- 


P 


(VOl 


(34) 


where  =  2T/H  or  the  ratio  of  the  Magnus -lift  coupling  to  the 

damping.  It  follows,  from  Equation  (33),  that  P/[<J>^  “ 

must  be  less  than  unity  for  a  statically  stable  missile  and  greater 
than  unity  lor  a  statically  unstable  missile.  In  either  case. 
Equation  (34)  may  be  rewritten  using  Equation  (33)  to  give: 


Conditions  for  stability  (A^,  A2  <  0)  are  that 


1 


> 


Vi. 

JF7s9 


or  equivalently 


> 


1  -  '/s, 


(35) 


The  boundary  between  stability  and  instability  follows  as 


(36) 


which  is  Equation  7.7  of  Reference  (2). 
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A  dynamically  stable  model  must  have  an  inverse  gyroscopic 
stability  parameter  less  than  that  given  by  Equation  (36) . 

If  a  trajectory  is  determined  for  a  given  projectile  and  if 
1/E-g  and  are  calculated  at  a  fixed  interval,  then  the  locus 

of  the  curve  in  the  1/Sg  -  plane  readily  indicates  the  stability 
over  its  flight  path. 

It  is  possible  now  to  write  1/Sg  and  in  terms  of  the 

aerodynamic  coefficients  and  projectile  mass,  mass  distribution  and 
geometry.  In  terms  of  problem  variables,  Sg  is: 

1  =  =  _ 2isi-fe - - —  (37  j 

or  in  terms  of  computer  variables. 


similarly, 


2.7  _ 

and  in  terms  of  computer  variables, 

D-2.*JS-.5*CKfe/XZ)*-C 


(37b) 


(38a) 


(38b) 


It  has  been  shown  that  the  solution  of  Equation  (30)  will  be 
stable  if  1/Sg  has  a  value  no  greater  than  that  given  by  Equation 
(36). 

The  third  term  on  the  right  in  Equation  (30)  is  the  steady- 
state  value  of  £,  this  term  is  designated  as  £g  and  called  the 
yaw  of  repose.  This  term  is,  according  to  Murphy  (Ref.  (2),  Eq.  6.13), 
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-P 


M  +•  \ 


(39) 


Inserting  the  appropriate  values  for  M,  P  and  T  from  Equations  (29) , 
the  above  expression  becomes,  after  rationalizing  the  denominators 


"X  9  CoS  & 


(40a) 


Since  pd/Vo  is  the  order  of  0.1,  the  imaginary  terms  in  Equation  (40a) 

can  be  ignored.  If  this  is  done,  it  is  possible  to  rewrite 
Equation  (40a)  in  terms  of  the  computer  variables  as, 

S3  _  (II*  PS*C,4.34fr*ftAA  (40b) 

(R  *bS*-CV**)*E>) 


In  calculating  Sg  (Eqs.  (37))  it  is  necessary  to  insert  a 
value  for  the  density.  The  density,  of  course,  changes  with 
altitude  and  this  change  is  assumed  in  this  report  to  be 
exponential  with  altitude,  y,  as: 


(41a) 


The  main  assumption  in  deriving  Equation  (41a)  is  that  the 
atmosphere  is  isothermal.  In  computer  variables,  Equation  (41a) 
becomes 


"R-  .O0Z377*  EXP(-.  3ife-?ZE-4  n't) 


(41b) 


The  aerodynamic  coefficients  which  are  necessary  to 
calculate  S^,  Sg  and  £g  are  enuered  into  the  program  as  a 

table  with  Mach  number.  Mach  number  is,  of  course,  the  ratio 
of  the  body's  speed  to  the  speed  of  sound  as 
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H=  V/C 


The  body's  speed  is  obtained  in  the  particle  trajectory  program, 
Equations  (17) .  The  speed  of  sound  varies  solely  with  the 
temperature  for  an  ideal  gas  and  in  the  English  units  may  be 
given  as, 


C  «■  ^3.0  h  +■  4  k° 


where  T  is  in  degrees  Fahrenheit  and  c  is  in  feet  per  second. 
The  speed  of  sound  may  be  written  as 


C  -V-  —  Vi  -  C  +  ii:  v i 

°  3*0  o  0 


where 


M  ■=  ^  /2C. 

aT 

from  Equation  (43)  and  3T/3y  (partial  derivative)  will  be  taken 
at  -0.0039  degrees  Fahrenheit  per  foot.  Thus,  Equation  (44)  may 
be  written  as 


c  «  i\ i 7  ->ooizy 


The  Mach  number  may  be  written  in  computer  variables  as 


V/CIU7  -  og4Z.*\) 


Equation  (44)  and,  hence,  Equation  (46)  are  assumed  to  be  valid 
only  up  to  36,000  feet.  Above  this  altitude  the  temperature 
and,  therefore,  the  speed  of  sound  are  assumed  to  be  constant. 
In  this  case,  Equation  (46)  will  be  replaced  by 

M “  v/q~M 
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Table  2  gives  the  correspondence  between  the  computer  variables , 
as  used  in  the  BASIC  language ,  and  the  problem  variables, 
most  of  which  have  been  defined  in  the  list  of  symbols.  The 
general  rules  of  the  BASIC  language  and  terminal  apply.  One  need 
only  supply  the  input  data  starting  on  line  1890  (see  Table  2) . 

Data  are  input  in  the  following  order: 

T5  initial  time  (seconds) 

XO  initial  horizontal  displacement  (feet) 

YO  initial  vertical  displacement  (feet) 

VO  platform  velocity  relative  to  inertial  space 

(feet  per  second) 

VI  projectile  muzzle  velocity  (feet  per  second) 

<f>  initial  flight  path  angle  (degrees) 

H  program  parameter  to  control  integration  step 

size  (seconds) 

N  program  parameter  to  control  the  print-out  frequency 

(integer  multiple  of  H) 

W  mass  of  projectile  (slugs) 

D5  diameter  of  projectile  (feet) 

11  axial  moment  of  inertia  (slug-ft2) 

12  transverse  moment  of  inertia  (slug-ft2) 

P5  initial  spin  rate  (radians  per  second) 

Z9  an  integer  designating  the  table  length  (number 

of  Mach  number  entries  -  maximum  value  is  ten) 

Following  the  above  data  are  the  tables.  The  tables  are 
supplied  per  line: 


M(J) 
A  (J) 

B  (J) 
C(J) 
D  (J) 


Mach  number 


C 


D 


C  + 
m  m 

q 


E(J)  C£ 


F  (J) 
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The  input 

data  should  appear  as 

follows: 

1890  DATA 

T5,  XO,  YO,  VO,  VI,  0 

1900  DATA 

H,  N 

1910  DATA 

W,  D5 ,  11,  12,  P5 

1920  DATA 

Z9 

1930  DATA 

M(l)  ,  A (1)  ,  B  (1)  , 

C(l)  ,  D  (1)  ,  E  (1)  ,  F  (1) 

1940  DATA 

M(2)  ,  A(2)  ,  B(2), 

C  (2)  ,  D  (2)  ,  E  (2)  ,  F  (2) 

19..  DATA  M(Z9)  ,  A(Z9),  B(Z9),  C(Z9),  C(Z9),  E(Z9) 


Typical  output  from  the  program  is  shown  in  Table  4.  As 
indicated,  the  first  part  of  the  output  is  the  coefficient  table. 

After  this  table  is  printed,  the  quantitites  H  and  N  are  printed. 

These  numbers  indicate  the  step  size  and  print-out  frequency, 
respectively.  The  remainder  of  the  print  out  is  the  trajectory 
parameters  and  stability  derivatives  (as  defined  in  Table  1) . 

These  data  are  presented  in  pairs.  Each  pair  contains  TIME 
(seconds) ,  X  (horizontal  distance  -  feet) ,  Y  (vertical  distance  - 
feet) ,  V  (speed  relative  to  inertial  space  -  feet  per  second) , 
p  (spin  rate  -  radians  per  second) ,  1/Sg  (inverse  gyroscopic 
stability  parameter) ,  Sd  (damping  stability  parameter) ,  Yaw 

of  repose  (radians) ,  pd/2V  (reduced  frequency) ,  MACH  (Mach  number) . 
Following  the  initial  conditions  is  a  print  out  of  the  above 
trajectory  constants  at  the  appropriate  step  dictation  by  the 
quantities  H  and  N.  The  program  automatically  terminates  at  zero 
altitude.  It  will  be  noted,  in  Table  3,  that  the  projectile 
velocity  relative  to  inertial  space,  V,  is  calculated  in  line 
430  from  the  platform  velocity  and  muzzle  velocity. 

STABILITY  ANALYSIS 

In  considering  the  stability  of  the  projectile  the  analyst 
must,  in  some  way,  limit  the  infinite  number  of  initial  conditions 
in  altitude  and  velocity  that  the  projectile  theoretically  might 
encounter.  A  relatively  small  number  of  initial  conditions  can 
be  chosen  which  will  represent  typical  operating  conditions.  Since 
the  20mm  projectile  under  consideration  is  intended  for  general 
purpose,  one  must  consider  firings  from  high-speed  aircraft, 
helicopters  and  ground. 

In  order  to  represent  firings  or  launchings  -Prom  high-speed 
aircraft,  the  gyroscopic  stability  factor  was  calculated  at  the  muzzle  for 
launchings  from  aircraft  flying  at  300,  475  and  600  knots  and  at 
altitudes  of  0,  5000,  10,000  and  15,000  feet.  In  addition,  since 
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atmospheric  conditions  are  influential,  it  was  decided  to  carry 
out  identical  calculations  for  a  cold  day  (-65°F)  and  a  standard 
day  (-59° F).  Using  the  program  described  in  the  preceding  section, 
the  gyroscopic  stability  factor  was  rapidly  and  easily  determined 
for  these  24  various  launchers. 

Figure  50  presents  results  for  the  cold  day  and  Figure  51 
for  the  standard  day.  A  comparison  of  these  figures  indicates 
that  the  most  critical  firing  condition  of  those  studied  occurs  on 
a  cold  day  in  firing  from  an  aircraft  having  a  6Q0  know  airspeed 
at  Sea  level.  Clearly,  the  projectile's  stability  decreases  with 
an  increase  in  speed  and  decrease  in  altitude  of  the  launch  aircraft 
and  with  decreasing  ambient  temperature.  Both  Figures  50  and  51 
indicate  that  the  projectile  is  gyroscopically  stable  for  all 
conditions  studied. 

It  was  stated  in  Equation  (36)  that  the  necessary  and  sufficient 
condition  for  projectile  dynamic  stability  is  that 


1  / s j  <  i  o 

and 


(48a) 


(48b) 


It  can  be  shown  that  a  good  approximation  to  Equation  (38a) 
might  be  written  neglecting  the  normal-force,  ,  and  the  drag,  C^, 

a 

contributions 


T  C 
nrv  gf 


(49) 


Now  since  Cm  +  Cm>  will  always  be  negative,  it  can  be  seen  that 
mq  a 

the  sign  of  will  be  identical  to  the  sign  of  the  Magnus  derivative, 

C  .  It  will  be  noted,  in  Figure  49,  that  the  Magnus  derivative 
mpct 

changes  sign  from  positive  to  negative  in  going  from  supersonic  to 
subsonic  flight.  Thus,  in  light  of  Equations  (48),  the  projectile 
will  be  unstable  for  subsonic  flight.  Figures  35  and  49  also  point 
out  that  there  is  considerable  scatter  in  both  the  damping-in-pitch 
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derivative  and  the  Magnus  derivative.  There  might  be  some  doubt  as 
to  the  projectile  satisfying  the  condition  that  <  2.  The 

stability  conditions,  as  specified  by  1/Sg  and  S^,  can  now  be 

examined  for  a  number  of  representative  operational  conditions. 

Since  projectile  stability  becomes  more  severe  the  higher 
the  platform  velocity,  it  was  decided  to  examine  the  trajectory 
of  the  projectile  fired  from  an  aircraft  flying  at  an  airspeed  of 
600  knots  and  at  several  altitudes.  In  Figures  52  through  55, 
the  reciprocal  of  the  gyroscopic  stability  factor  is  plotted 
versus  the  damping  stability  factor  for  launchings  from  an 
aircraft  having  an  airspeed  of  600  knots  and  altitudes,  respectively, 
of  15,000,  10,000,  5,000  and  500  feet.  Certain  conclusions 
about  projectile  stability  may  now  be  drawn  from  these  figures. 

In  Figure  52  it  will  be  noted  that  the  projectile  remains 
gyroscopically  stable  throughout  the  entire  flight.  This  conclusion 
was  made  on  the  basis  of  Figures  50  and  51.  More  significantly, 
the  dynamic  stability  factor  decreases  steadily  with  diminishing 
flight  speed.  Thus,  an  error  made  in  the  location  of  due  to 

scatter  is  not  too  serious  as  the  damping  stability  factor  decreases 
steadily  through  the  flight.  In  comparing  Figures  52  through  55 
it  will  be  noted  that  stability  becomes  more  critical  at  the 
lower  altitudes,  not  only  gyroscopic  stability,  as  pointed  out 
earlier,  but  damping  stability  as  well.  It  is  interesting  to  note 
that  at  15,000  feet  the  unstable  conditions  are  encountered  at 
between  six  to  seven  seconds  after  firing  and  around  13,000  feet 
downrange.  By  contrast,  at  500  feet  instability  is  encountered 
at  less  than  four  seconds  and  at  a  downrange  distance  of  about 
7500  feet.  In  both  extremes,  projectile  velocities  are  comparable 
at  the  onset  of  instability. 

Since  the  projectile  under  consideration  is  general  purpose, 
is  was  decided  to  examine  its  stability  under  conditions  of 
horizontal  firing  from  a  stationary  platform  and  a  45-degree 
firing  from  the  ground.  From  the  stationary  platform  the  instability 
is  encountered  at  less  than  three  seconds  and  about  5,000  feet  down- 
range.  For  ground  firing  at  45  degrees,  the  instability  occurs 
at  an  altitude  and  range  both  of  about  3500  feet  (see  Figs.  (56) 
and  (57) ,  respectively) . 


CONCLUSIONS 

This  report  has  presented  aerodynamic  data  on  the  U.  S.  Navy's 
general  purpose  20-millimeter  projectile.  The  data  were  obtained 
in  both  a  pressurized  aeroballistics  range  and  in  a  supersonic 
wind  tunnel,  h  computer  analysis  was  carried  out  to  determine 
the  stability  of  the  projectile  under  several  operational 
situations.  It  may  be  concluded  that  the  projectile  is  both 
gyroscopically  and  dynamically  stable  under  all  conditions 
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examined.  Gyroscopic  stability  becomes  lessened  with  decreasing 
altitude  and  temperature  and  increasing  platform  velocity.  The 
dynamic  stability  decreases  with  decreasing  altitude,  although 
it  is  relatively  insensitive  to  platform  velocity. 

In  addition  to  the  above  computer  stability  analysis, 
attempts  were  made  to  calculate  the  normal-force  and  pitching- 
moment  derivatives,  and.C^  ,  respectively,  using  Wood's 

a  a 

semiempirical  method.  The  agreement  is  considered  only  fair. 
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TABLE  1 

PROGRAM  FOR  CALCULATION  OF  NORMAL  FORCE  AND  PITCHING  MOMENT  DERIVATIVES  FOR  A  SHELL 

100  HEM  THIS  PROGRAM  CALCULATES  NORMAL  FORCE  AND  HITCHING  MOMENT 
110  REM  DERIVATIVES  AND  CENTER  0F  PRESSURE  Fv3r  A  SHELL. 

120  REM  THIS  IS  A  SEMI -EMPIRICAL  METHOD  BASED  UPON  BKL854 
130  REM  Cl  IS  THE  NORMAL  FORCE  DERIVATIVE 

1  AO  REM  C2  IS  THE  PITCHING  MOMENT  DERIVATIVE  ABOUT  THE  BASE 

150  REM  C3  IS  THE  PITCHING  MOMENT  ABOUT  THE  CENTER  OF  GRAVITY 

160  REM  S  IS  THE  STATIC  MARGIN 

170  REM  CO  IS  THE  CENTER  OF  PRESSURE  IN  CALIBERS  FROM  THE  NOSE 
130  REM  R  IS  THE  FINENESS  RATIO  »  BODY  VOLUME/CrEF.  AKEAXDI AMETEr) 

190  REM  G  IS  THE  LOCATION  OF  THE  CG  FROM  THE  BASE 

200  REM  L  IS  THE  LENGTH  OF  THE  BODY  IN  CALIBERS 

210  REM  T  IS  THE  CENTER  OF  PRESSURE  AFT  OF  THE  VERTEX  (IN  CAL.) 

220  REM  A 1  IS  THE  RATIO  OF  EFFECTIVE  BOAT TAI L  CROSSECTI ONAL  AREA 

230  REM  TO  REFERENCE  AREA* (AB+AR)/2*AK  WHERE  AB  IS  BOATTAIL  AkEA 

240  REM  AND  AR  IS  THE  KEFFERENCE  AREA 

250  FOR  I  =  1  TO  15 

260  READ  M 

270  LET  R  =  3.83148 

280  LET  L  =  4.9 

290  LET  A 1  =1 

300  LET  N0=. 14331 

310  LET  Nl=-9.32777E-2 

320  LET  N2=. 453401 

330  LET  N3= .581451 

340  LET  N4=- 4 .97192 

350  LET  N5=6. 97863 

360  LET  N6=-3. 80723 

370  LET  N7=. 736042 

380  LET  M0=. 192947 

390  LET  M 1 =- • 1 12361 

400  LET  M2=. 872645 

410  LET' M3=-2. 37587 

420  LET  M4*l .331  14 

430  LET  M5=. 924154 

440  LET  M6=- 1.07323 

450  LET  M7=. 26358 

460  LET  G  =  2.2  1 

470  IF  M> 1  THEN  510 

480  LET  FI  =  l 

490  LET  F2  =  1 

500  GOTO  580 

510  LET  N  =  SQR((M»2)- 1)/R 

520  LET  Fl=N0+Nl*N+N2*CNT2)+N3*(N»3>+N4*CNt4>+NS*(Nt5) 

530  LET  F  1  =F 1 +N6* ( N  t  6 ) +N7* ( N 1 7 ) 

540  LET  F2=M0+M 1*N+M2* (N  t  2) +M3* (N »  3>  +  M4* ( N  f  4 >  +  M5* ( N t  5) 

550  LET  F2=F2+M6*(N»6)+M7*(Nt 7) 

560  LET  FI  a  F1*SQRC(Mt2)-1)+1 

570  LET  F2  =  (F2)*SQRC  <Mt2)- 1  )■*■  1 

580  LET  Cl  =  ( 2*A 1  +  »  5) *F 1 

590  LET  C2=2*R*F2 

600  LET  C=C2/C 1 

610  LET  T  =  L  -C 

620  LET  S=C-G 

630  LET  C3=C 1 *S 

640  IF  I>1  THEN  670 

650  PRINT  "PITCH.  MOM.  "*"N0R.  FORCE"*"STAr.  MAR.  "* "C . P. "* "MACH  NO." 

660  PRINT 

670  PRINT  C3*C 1  *  S*  T* M 
680  NEXT  I 

690  DATA  .4* i .* 1 .2* 1 . 5* 1 • 75*2* 2. 5* 3* 3. 5 
700  END 
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Table  2 


CORRESPONDENCE  BETWEEN  COMPUTER  AND  PROBLEM  VARIABLES 


Computer 

Problem 

Variable 

Variable 

Definition 

A 

CD 

drag  coefficient,  current  value 

A1 

D 

drag,  CDpSV2/2 

A(Z) 

— 

drag  coefficient  array 

B 

Cm 

pitching-moment  derivative,  current  value 

B  (2) 

cx 

pitching-moment  array 

C 

c  +  c 

m  m. 

q  a 

damping-in-pitch  derivative,  current  value 

C(Z) 

— 

damping-in-pitch  derivative  array 

D 

CN 

a 

normal-force  derivative 

D1 

variable  used  in  Runge  Kutta 

D2 

— 

variable  used  in  Runge  Kutta 

D5 

d 

diameter  of  projectile,  reference  length 
(feet) 

D(Z) 

— 

normal-force  derivative  array 

E 

C£ 

p 

damping-in-roll  derivative 

El 

variable  used  in  Runge  Kutta 

E2 

— 

variable  used  in  Runge  Kutta 

E5 

— 

exponential  part  of  Pd/2V  formula 

E(Z) 

-- 

damping-in-roll  derivative  array 

F 

c 

n  „ 
pa 

Magnus  moment  derivative 

F(Z) 

— 

Magnus  moment  derivative  array 

G 

— 

counter 

G5 

— 

switching  variable 

H 

-- 

step  size 

I 

— 

counter 

11 

I 

XX 

2 

axial  moment  of  inertia  (slug-ft  ) 

12 

I 

yy 

2 

transverse  moment  of  inertia  (slug-ft  ) 
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Table  2  (Cont'd) 


Computer 

Problem 

Variable 

Variable 

Definition 

J 

— 

index  variable 

K1 

variable  used  in  Runge  Kutta 

K2 

— 

variable  used  in  Runge  Kutta 

K3 

— 

variable  used  in  Runge  Kutta 

K4 

— 

variable  used  in  Runge  Kutta 

K6 

md2 

mass  times  reference  length  squared 
(slug-ft2) 

LI 

— 

variable  used  in  Runge  Kutta 

L2 

— 

variable  used  in  Runge  Kutta 

L3 

— 

variable  used  in  Runge  Kutta 

L4 

— 

variable  used  in  Runge  Kutta 

M 

— 

Mach  number 

Ml 

— 

variable  used  in  Runge  Kutta 

M2 

— 

variable  used  in  Runge  Kutta 

M3 

— 

variable  used  in  Runge  Kutta 

M4 

— 

variable  used  in  Runge  Kutta 

M(Z) 

— 

Mach  number  array 

N 

— 

print-out  frequency 

N1 

— 

variable  used  in  Runge  Kutta 

N2 

— 

variable  used  in  Runge  Kutta 

N3 

— 

variable  used  in  Runge  Kutta 

N4 

— 

variable  used  in  Runge  Kutta 

0 

*o 

initial  flight  path  angle  in  degrees 

P 

y 

dy/dt  velocity  in  vertical  direction 

y 

in  ft/sec 

PO 

— 

variable  used  in  Runge  Kutta 

PI 

— 

dP/dt,  P,  dVy/dt 

P5 

p 

current  value  of  spin  rate  in  rad/sec 

P7 

— 

nonexponential  part  of  pd/2V  formula 

P8 

pd/2V 

current  value  of  pd/2V 

Q 

Vx 

dx/dt,  velocity  in  horizontal  direction 
in  ft/sec 
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Computer 

Variable 

QO 

Ql 

Q5 

R 

S 

51 

52 

53 
T 

TO 

T5 

V 

VO 

VI 

V2 

V5 

W 

X 

XO 

XI 
X5 

Y 

YO 

Yl 

Z 

Z9 


Table  2  (Cont'd) 

Problem 

Variable  Definition 

variable  used  in  Runge  Kutta 
dQ/dt,  Q,  dVk/dt 

Vxo  previous  value  of  Vx  in  ft/sec 

p  density  slugs/ft2 

S  reference  area  Trd2/4  in  ft2 

Sg  gyroscopic  stability  factor 

S^  dynamic  stability  factor 

6y  yaw  angle  of  repose  in  rad 

t  time  in  seconds 

variable  used  in  Runge  Kutta 

to  initial  time  in  seconds 

v  velocity  of  projectile  relative  to 

inertial  space 

velocity  of  platform  relative  to 
inertial  space 

muzzle  velocity 

v2  total  velocity  squared  in  ft^/sec 

vo  previous  value  of  V  in  ft/sec 

m  mass  of  projectile  in  slugs 

x  horizontal  distance  in  feet 

variable  used  in  Runge  Kutta 
dx/dt  horizontal  velocity 

Xo  previous  value  of  X  in  feet 

y  vertical  distance  in  feet 

variable  used  in  Runge  Kutta 
dy/dt  vertical  velocity 

index  variable 
table  length 
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TABLE  3 

BASIC  PROGRAM  FOR  THE  CALCULATION  OF  PROJECTILE  STABILITY 

500  REM  THE  PURPOSE  OF  THIS  PROGRAM  IS  TO  OBTAIN  THE  STABILITY 
110  REM  OF  A  PROJECTILE  IN  TERMS  OF  THE  RECIPROCAL  OF  THE 
120  REM  GYROSCOPIC  STABILITY  FACTOR  AND  THE  DAMPING  STABILITY 
130  REM  FACTOR.  THE  FOLLOWING  IS  A  LISTING  OF  THE  REQUIRED 
140  REM  INPUTS  TO  THIS  PROGRAM. 

ISO  REM  T5  IS  THE  INITIAL  TIME(SECONDS) 

160  REM  XO  IS  THE  INITIAL  HORIZONTAL  DI SPLACEMENTCFEET) 

170  REM  YO  IS  THE  INITIAL  VERTICAL  DI SPLACEMENTCFEET) 

180  REM  V  IS  THE  VEL.  OF  PKOJ.  KEL.  TO  INERTIAL  SPACEC F T/SEC . ) 

190  REM  VO  IS  THE  PLATFORM  VELOCITY  rEl.  TO  INERT.  SPACE  CFT/SEC) 

200  REM  VI  IS  THE  PROJECTILE  MUZZLE  VELOCITY  ) FT/SEC) 

210  REM  0  IS  THE  INITIAL  FLIGHT  PATH  ANGLEC DEGREES) 

220  REM  H  IS  PROGRAM  PARAMETER  TO  CONTROL  I  NT.  STEP  SIZECSEC.) 

230  REM  N  IS  PROGRAM  PARAMETER  TO  CONTROL  PRINT-OUT  FREQUENCY 

240  REM  (INTEGER  MULTIPLE  OF  H) 

250  REM  W  IS  THE  MASS  OF  THE  PROJECTILE  (SLUGS) 

260  REM  D5  IS  THE  DIAMETER  OF  THE  PROJECTI LE(FEET) 

270  REM  II  IS  THE  AXIAL  MOMENT  OF  I NERTI AC SLUG-FT  SQ.) 

280  REM  12  IS  THE  TRANSVERSE  MOMENT  OF  INERTIA  (SLUG-FT  SQ.) 

290  REM  Z9  IS  AN  INTEGER  INDICATING  THE  LENGTH  OF  THE  TABLES 

300  REM  (NUMBER  OF  MACH  NUMBER  ENTRI ES-MAXI MUM  NO.  IS  10) 

310  REM  THE  TABLES  ARE  SUPPLIED  PER  LINE  AS  FOLLOWS: 

320  REM  M(J)  IS  THE  MACH  NUMBER 

330  REM  A(J)  IS  THE  DRAG  COEFFI Cl  ENT, CD 

340  REM  BCJ)  IS  THE  PITCH.  MOMENT  DERIVATIVE*  CMA 

350  REM  CCJ)  IS  THE  DAMPING- IN- PITCH  DERI VATI VE> CMG+CMA 

360  REM  D(J)  IS  THE  NORMAL  FORCE  DERI VATI VE, CNA 

370  REM  E(J)  IS  THE  DAMPING-IN-ROLL  DERI VATI VE, CLP 

380  REM  F ( J )  IS  THE  MAGNUS  DERI VATI VE, CMPA 

390  PRINT  »  PROJECTILE  STABILLTY  PROGRAM" 

400  PRINT 

410  READ  T5»X0,Y0, VO, V1,0,H,N 
420  LET  V  =  VO  +  VI 
430  READ  W,D5, I 1,I2,P5 

440  DIM  MC10),  AC  10),  8(10),  C(10),  D(10),EC10),  F(10> 

450  READ  Z9 

460  PRINT  "COEFFICIENT  TABLE" 

470  PRINT  "  M  A  B  C  D  E  F" 

480  FOR  J=1  TO  Z9 

490  READ  M ( J ) ,  A(J),  BCJ),  C(J),  D(J),  E(J)»  FCJ) 

500  PRINT  M(J) )A(J)sB(J))C(J> I D( J) J E( J) 8 F( J) 

510  NEXT  J 
520  LET  Z=1 

530  LET  S=(3.14159*D5»2)/4 
540  LET  P0=V*S IN (0/57.296) 

550  LET  QO=V*C9S(0/57 .296) 

560  LET  K 6  = ( W* ( D 5 * 2 ) ) 

570  LET  T0=T5 
580  PRINT 

590  PRINT  "FOR  H  ="HJ"N  ="N 
600  PRINT 
610  PRINT 

620  PRINT  "TIME", "X", "Y"> "V", "P" 

630  PRINT 

640  PRINT  " 1 /SG"» "SD", "YAW  OF  REPOSE", "PD/2V”, "MACH" 

650  PRINT 
660  PRINT 

670  PRINT  "INTITIAL  CONDI  TONS: " 

630  PRINT 
690  LET  T=TO 
700  LET  T0=T5 
710  LET  X=XO 
720  LET  Y=YO 
730  LET  G=QO 
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(CONT'D) 


740  LET  P=P0 

750  LET  X5=X 

760  LET  US=U 

770  LET  V5=V 

780  GOSUB  1620 

790  LET  G5=0 

800  LET  G=0 

810  LET  I  =  N 

820  LET  G=G+  1 

830  IP  G5=0  THEN  MOO 

8-10  LET  G5=  1 

850  GOSUB  1560 

860  LET  Kl=Yl*H 

870  LET  N 1 =X 1 *H 

880  GOSUB  1590 

890  LET  L1=P1*H 

900  LET  M1=01*H 

910  LET  T  =T0+H/2 

920  LET  Y  =Y 0+K 1/2 

930  LET  X=XO+N 1/2 

940  LET  P=PO+L 1/2 

950  LET  G=C0+Ml/2 

960  GOSUB  1560 

970  LET  K2=Y1*H 

980  LET  N2=X 1 *H 

990  GOSUB  1590 

1000  LET  L2=P1*H 

1010  LET  M2=QI*H 

1020  LET  Y=Y0+K2/2 

1030  LET  X=X0+N2/2 

1040  LET  P=P0+L2/2 

1050  LET  Q=Q0+M2/2 

1060  GOSUB  1560 

1070  LET  K3=Y1*H 

1080  LET  N3=X 1 *H 

1090  GOSUB  1590 

1100  LET  L3=P1*H 

1110  LET  M3=0 1 *H 

1120  LET  T=TO+H 

1 130  LET  Y= Y0+K3 

1140  LET  X=X0+N3 

1150  LET  P=  P0+L3 

1160  LET  0=U0+M3 

1 170  GOSUB  1560 

1180  LET  K4=Y 1*H 

1190  LET  N4=X1*H 

1200  GOSUB  1590 

1210  LET  L4=P1*H 

1220  LET  M4=0 1 *H 

1230  LET  Dl=(Kl+K4)/6+<K2+K3)/3 

1240  LET  El=(Nl+N4)/6+<N2+N3)/3 

1250  LET  D2=(Ll+L4)/6+(L2+L3)/3 

1260  LET  E2=(Ml+M4)/6+(M2+M3)/3 

1270  LET  Y0= YO+D 1 

1280  LET  X0=X0+E 1 

1250  LET  P0=P0+D2 

1300  LET  00=Q0+E2 

1310  LET  T0=T5+H*G 

1320  LET  T=TO 

1330  LET  Y=Y0 

1340  LET  X=XO 

1350  LET  P=PO 

1360  LET  0=00 

1370  GOSUB  1620 
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(CONT'D) 

1380  IF  Y<  =  0  THEN  1400 
1390  IF  G<I  THEN  820 

1400  LET  B=BCZ>-KCB(Z+l)-B(Z>>*(M-M<Z5>>/CMCZ+!>-MCZ>> 

1410  LET  C=GCZ>«-(CC<Z+1)-C(Z>>*(M-M(Z>>>/(M<Z+1)-MCZ>> 

1420  LET  D=D<Z>+CCD<Z+1>-DCZ>)*<M-M(Z>>>/CMCZ+1>-MCZ>) 

1430  LET  F=FCZ>+{(F<Z+1>-F(Z))*<M-MCZ>>>/<M<Z+1>-M<Z>> 

1440  LET  Sl=C2*(Ilf2)*(P5t2))/(I2*3»1416*K*(D5»3>*(V»2)*B> 

1450  LET  SI=1/S1 

1460  LET  S2=<2*(D-A+.5-«<K6/I  1 )  *F>  >  / ( D- <  2* A)  -  .  5* < K 6/1 2>  *C > 

1470  LET  S3=(  I  1*P5*64  «348*C 0/V))/(R*S*D5*CVt3)*B> 

1480  PRINT  T,X,Y,V,P5 
1490  PRINT  S 1 ,  S2» S3,  P8,  M 
1500  PRINT 
1510  PRINT 

1520  IF  Y<=0  THEN  2030 
1530  IF  G5=0  THEN  840 
1540  LET  I  =  I  +  N 
1550  GOTO  820 
1560  LET  X1=0 
1570  LET  Y 1 =P 
1580  RETURN 

1590  LET  Pl=-((Al)/(vJ*V))*P-32.  174 

1600  LET  Ul=-C  (A1  >*0 

1610  RETURN 

1620  LET  V2=P*P+U*Q 

1630  LET  V=S0R(V2) 

1640  IF  Y <3 6 500  THEN  1670 

1650  LET  M=V/971 

1660  GO  TO  1680 

1670  LET  N=V/( 1117-.0042+Y) 

1680  IF  M<=M(Z+ 1 )  THEN  1720 

1690  LET  Z=Z+l 

1700  IF  Z>  Z9  THEN  1880 

1710  GO  TO  1680 

1720  IF  M>=M(Z)  THEN  1760 

1730  LET  Z=Z-1 

1740  IF  Z<1  THEN  1880 

1750  GO  TO  1720 

1760  LET  A=A(Z >  +  ((A(Z+1)-A(Z)) * (M-M(Z) >)/<M<Z+l)-M<Z)) 

1770  LET  H=  •  0023 7 7* EXP ( - • 3 1 582E-4* Y> 

1780  LET  Al=(A*S*R*V2>/2 

1790  LET  E=E(Z)  +  ((E(Z+1)-E<Z> ) *<M-M(Z) ) )/(M( Z+ 1 ) -N< Z) ) 

1800  LET  P7=<P5*V)/V5 

1810  LET  E5=£XP( { V5/Q5)*<  <  (R*S)/(2*W)  )*{  (K6/I  1 >* . 5*E+A)*(X-X5) > > 

1820  LET  P5=P7*E5 

1830  LET  P8=(P5*D5)/(2*V) 

1840  LET  X5=X 
1850  LET  05=0 
I860  LET  V5=V 
1870  RETURN 

1880  PRINT  "ERROR  MACH  NO.  =  "M"  NOT  IN  DATA  " 

1890  DATA  0,0,1000,1010,2700,0 
1900  DATA  .1.10 

1910  DATA  .0100578. .06562. .5283E-5. .601331E-4. 10285 
1920  DATA  10 

1930  DATA  0. .25.2.5.-22.5. 1.-. 044.-1 
1940  DATA  .8. .275. 2. 84, -22. 5. 2. 19. -.039. -.76 
1950  DATA  1, .435, 2. 98, -22. 5, 2. 42, -.0375,-. 45 
1960  DATA  1 . 1, .515,3. 06, -22. 5, 2. 52,-. 037, -.2 
1970  DATA  1.2, .523, 3. 09, -22. 5, 2. 63, -.0365,. 15 
1980  DATA  1 .5, .5,3. 13, -22. 5, 2. 9, -.0345, .75 
1990  DATA  2, .412,2.92, -22. 5,3. 22,-. 0315, 1 . 15 
2000  DATA  2.5, .33, 2. 48, -22. 5, 3. 5, -.0285, 1.2 
2010  DATA  3, .285, 1 .88, -22. 5,3.4, - .025, 1 .3 
2020  DATA  5, .275, 1 .88, -22 . 5, 3 . 4, - . 025, 1 .35 
2030END 
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TABLE  -'i  SAMPLE  OUTPUT  FROM  BASIC  STABILITY  PROGRAM 

PKJJEC  TILL  STABILITY  PKOGkAM 


OJKKrIClfcAlf  f  AoLfc. 


•1  A 

a  C 

0  E 

F 

0  2.5  - 

1  -.04/4 

-  1 

•A  -2/5  2*88  -22.5  2.19 

-.039  -.76 

1  • 435  2.9B 

-22.5  2.82  - 

.0375  -.85 

1  •  1  .bib  1 . 

06  -22.5  2.52 

-.037  -.2 

1.2  -523  3. 

09  -22.5  2.63 

-.0365  .  15 

1.5  .5  3.13 

-22.5  2.9  -. 

0385  .75 

2  .8  13  9.92 

-22.5  3.22  - 

.0315  1.15 

2.5  .33  2.88  -22.5  3.5  - 

•0285  1.2 

1  .265  l.HH 

-22.5  3.8  -. 

025  1.3 

5  .275  1.88 

-22.5  3.8  -. 

025  1.35 

FJr!  11=  .1  M  = 

10 

TIME 

Y 

9 

P 

l/SG 

.50 

YA.v  OF  KEPJSE 

PU/2V 

rtACH 

i.Mi  in  al  c.).9u  i  r:i9S: 

0 

U 

1000 

37  10 

10285 

.538756 

1 .55077 

7. 12572E-5 

9.09  57  IE-2 

3.33393 

1  . 

3086-87 

985.713 

2560.05 

9033.62 

. 869929 

1 .8607 

1 .38782E-8 

. 1 15776 

2.30083 

2. 

5280-29 

988.669 

1800.31 

8129.82 

.333251 

1  . 17706 

3 . 00232E- 8 

.1831 63 

1.6175 

3. 

6761.79 

898.587 

1326.  18 

7885.86 

•216828 

.538758 

6.B8708E-8 

. 188213 

1 .19127 

8 . 

7958.8 1 

828.715 

1063.08 

6898.39 

.  1  5508.-! 

-  .033608 

1  .28888E-3 

.212915 

.95865 

5. 

8988. 59 

736.226 

929.799 

6831  . 

. 1 30°26 

-.  1989  19 

1  .3  862E-3 

.228696 

.8387  18 

6  . 

9823.76 

626. 786 

886.668 

6150.03 

.  1  1  9228 

- .270187 

2.33566E-3 

.233327 

.759771 

7  . 

10625.8 

89  6.03  I 

760.785 

5875.88 

. 1 1058 

- .305803 

2.88106E-3 

.286913 

.700308 

8  . 

1  1363.2 

385.051 

''26.828 

5688 . 65 

.  1038 

-.33578! 

3.369 13E-3 

.258989 

•65118 

9  . 

12085.2 

1 78.379 

68 1 .208 

5951  • 

7.78I68E-2 

- .361018 

3 .905 1 9  E- 3 

.262585 

.610255 

10. 

12679.2 

-  15.2999 

683.823 

5290.76 

7.231S5E-J 

-.382819 

8.83269E-3 

.26979  1 

.575995 
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MAGNETIC  BALLBEARING 

TACHOMETER  \ 


MAGNUS  FLEXURE  MULTI-STAGE 

AIR  TURBINE 


FIG.  4  MAGNUS  BALANCE  DETAILS 
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GAGE  IN  COMPRESSION 


FIG.  6  MAGNUS  BRIDGE  UNDER  YAW  MOMENT 
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NORMAL  FORCE  COEFFICIENT 
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FIG.  10  NORMAL  FORCE  AND  PITCHING  MOMENT  COEFFICIENTS  VERSUS  ANGLE  OF  ATTACK  AT 
A  MACH  NUMBER  OF  0.79  FOR  THE  BASIC  CONFIGURATION 


FIG.  11  NORMAL  FORCE  AND  PITCHING  MOMENT  COEFFICIENTS  VERSUS  ANGLE  OF  ATTACK  AT 
A  MACH  NUMBER  OF  0.90  FOR  THE  BASIC  CONFIGURATION 


r.  s  j  tn  :t  I  Hi  t  H  » il  ifii  I  t !  !  J  I  i  { ft 1  il  1  i 


iniuiHuititatiaaioji 


FIG.  12  NORMAL  FORCE  AND  PITCHING  MOMENT  COEFFICIENTS  VERSUS  ANGLE  OF  ATTACK  AT 
A  MACH  NUMBER  OF  0.93  FOR  THE  BASIC  CONFIGURATION 
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FIG.  14  NORMAL  FORCE  AND  PITCHING  MOMENT  COEFFICIENTS  VERSUS  ANGLE  OF  ATTACK 
A  MACH  NUMBER  OF  1 .76  FOR  THE  BASIC  CONFIGURATION 


FIG.  15  NORMAL  FORCE  AND  PITCHING  MOMENT  COEFFICIENTS  VERSUS  ANGLE  OF  ATTACK 
A  MACH  NUMBER  OF  2.0  FOR  THE  BASIC  CONFIGURATION 
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FIG.  16  NORMAL  FORCE  AND  PITCHING  MOMENT  COEFFICIENTS  VERSUS  ANGLE  OF  ATTACK  AT 
A  MACH  NUMBER  OF  2.5  FOR  THE  BASIC  CONFIGURATION 
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FIG.  17  NORMAL  FORCE  AND  PITCHING  MOMENT  COEFFICIENTS  VERSUS  ANGLE  OF  ATTACK  AT 
A  MACH  NUMBER  OF  3.0  FOR  THE  BASIC  CONFIGURATION 
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FIG.  18  NORMAL  FORCE  AND  PITCHING  MOMENT  COEFFICIENT  VERSUS  ANGLE  OF  ATTACK  AT 
A  MACH  NUMBER  OF  0.80.FOR  THE  SMOOTH  CONFIGURATION 
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FIG  19  NORMAL  FORCE  AND  PITCHING  MOMENT  COEFFICIENTS  VERS! 
A  MACH  NUMBER  OF  0.95  FOR  THE  SMOOTH  CONFIGURATION 
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FIG.  20  NORMAL  FORCE  AND  PITCHING  MOMENT  COEFFICIENTS  VERSUS  ANGLE  OF  ATTACK  AT 
A  MACH  NUMBER  OF  1 .05  FOR  THE  SMOOTH  CONFIGURATION 


FIG.  21  NORMAL  FORCE  AND  PITCHING  MOMENT  COEFFICIENTS  VERSUS  ANGLE  OF  ATTACK  AT 
A  MACH  NUMBER  OF  1 .76  FOR  THE  SMOOTH  CONFIGURATION 

(p! 
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FIG.  22  NORMAL  FORCE  AND  PITCHING  MOMENT  COEFFICIENTS  VERSUS  ANGLE  OF  ATTACK  AT 
A  MACH  NUMBER  OF  2.0  FOR  THE  SMOOTH  CONFIGURATION 


FIG.  23  NORMAL  FORCE  AND  PITCHING  MOMENT  COEFFICIENTS  VERSUS  ANGLE  OF  ATTACK  AT 
A  MACH  NUMBER  OF  2.5  FOR  THE  SMOOTH  CONFIGURATION 

(P'b 


FIG.  24  NORMAL  FORCE  AND  PITCHING  MOMENT  COEFFICIENTS  VERSUS  ANGLE  OF  ATTACK  AT 
A  MACH  NUM8ER  OF  3.00  FOR  THE  SMOOTH  CONFIGURATION 
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FIG.  26  PIT 


4.0- 
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CENTER  OF  GRAVITY 
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FIG.  29  CENTER  OF  PRESSURE  FROM  BODY  VERTEX  VERSUS  MACH  NUMBER  FROM  BALLISTIC  RANGE 
AND  WIND  TUNNEL  MEASUREMENTS  FOR  THE  BASIC  CONFIGURATION 


NORMAL  FORGE  DERIVATIVE  VERSUS  MACH. NUMBER  FROM  BALLISTIC  RANGE  AND  WIND  TUNNEL 
MEASUREMENTS  FOR  THE  SMOOTH  CONFIGURATION 
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FIG.  3!  PITCHING  MOMENT  DERIVATIVE  VERSUS  MACH  NUMBER  FROM. BALLISTIC  RANGE 
AND  WIND  TUNNEL  MEASUREMENTS  FOR  THE  SMOOTH  CONFIGURATION 
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MEASUREMENTS  FOR  THE  BASIC  CONFIGURATI 
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MAGNUS  MOMENT  COEFFICIENT,  C  ;  MAGNUS  FORCE  COEFFICIENT 
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FIG.  38  MAGNUS  FORCE  AND  MOMENT  COEFFICIENTS  VERSUS  ANGLE  OF  ATTACK 
AT  A  MACH  NUMBER  OF  1.76  FOR  THE  BASIC  CONFIGURATION 
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FIG.  39  MAGNUS  FORCE  AND  MOMENT  COEFFICIENTS  VERSUS  ANGLE  OF  ATTACK 
AT  A  MACH  NUMBER  OF  2.0  FOR  THE  BASIC  CONFIGURATION 


MAGNUS  MOMENT  COEFFICIENT.  C  ;  MAGNUS  FORCE  COEFFICIENT 


ANGLE  OF  ATTACK  (DEGREES) 
O  C 


□  C 


42  MAGNUS  FORCE  AND  MOMENT  COEFFICIENTS  VERSUS  ANGLE  OF  ATTACK 
AT  A  MACH  NUMBER  OF  0.80  FOR  THE  SMOOTH  CONFIGURATION 


MAGNUS  MOMENT  COEFFICIENTS,  C  ;  MAGNUS  FORCE  COEFFICIENTS 
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FIG.  43  MAGNUS  FORCE  AND  MOMENT  COEFFICIENTS  VERSUS  ANGLE  OF  ATTACK 
AT  A  MACH  NUMBER  OF  0.93  FOR  THE  SMOOTH  CONFIGURATION 
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AT  A  MACH  NUMBER  OF  1.76  FOR  THE  SMOOTH  CONFIGURATION 


( 


AT  A  MACH  NUMBER  OF  2.0  FOR  THE  SMOOTH  CONFIGURATION 
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FIG.  47  MAGNUS  FORCE  AND  MOMENT  COEFFICIENTS  VERSUS  ANGLE  OF  ATTACK 
AT  A  MACH  NUMBER  OF  3.0  FOR  THE  SMOOTH  CONFIGURATION 


A  WIND  TUNNEL 


FIG.  49  MAGNUS  MOMENT  DERIVATIVE  VERSUS  MACH  NUMBER  FROM  BALLISTIC  RANGE 
AND  V/IND  TUNNEL  MEASUREMENTS  FOR  THE  BASIC  CONFIGURATION 


2„4- 


NOLTR  71-95 


AIRCRAFT  SPEEDS  OF  300,  475  AND  600  KNOTS 


FIG.  5!  GYROSCOPIC  STABILITY  FACTOR  VERSUS  ALTITUDE  FOR  AIR  TEMPERATURE  OF  59°  FARENHEIT  AND 
AIRCRAFT  SPEEDS  OF  300,  475,  AND  600  KNOTS 
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FIG.  52  l/Sg  VERSUS  Sd  FOR  THE  BASIC  BODY  LAUNCHED  HORIZONTALLY  FROM  AN  AIRCRAFT  FLYING  AT 
600  KNOTS  AIRSPEED  AND  15,000  FEET  ALTITUDE 


FIG.  53  1/Sg  VERSUS  S  ,  FOR  THE  BASIC  BOCY  LAUNCHED  HORIZONTALLY  FROM  AN  AIRCRAFT  FLYING 

a 

AT  600  KNOTS  AIRSPEED  AND  10,000  FEET  ALTITUDE 


FIG.  54  1/Sg  VERSUS  Sd  FOR  THE  BASIC  BODY  LAUNCHED  HORIZONTALLY  FROM  AN  AIRCRAFT  FLYING 
AT  600  KNOTS  AIRSPEED  AND  5000  FEET  ALTITUDE 


FIG.  55  1/Sg  VERSUS  Sd  FOR  THE  BASIC  BODY  LAUNCHED  HORIZONTALLY  FROM  AN  AIRCRAFT  FLYING 
AT  AN  AIRSPEED  OF  600  KNOTS  AND  500  FEET  ALTITUDE 
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FIG.  57  1/Sg  VERSUS  S  ,  FOR  BASIC  BODY  LAUNCHED  AT  GROUND  LEVEL  AT  45  DEGREE  ANGLE 


